Constraints on massive gravity theory from big bang nucleosynthesis 
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The massive gravity cosmology is studied in the scenario of big bang nucleosynthesis. By making 
use of current bounds on the deviation from the fractional mass, we derive the constraints on the 
free parameters of the theory. The cosmological consequences of the model are also discussed in the 
framework of the PAMELA experiment. 

PACS numbers: 



I. INTRODUCTION 



Recently, de Rham, Gabadadze and Tolley (dRGT) [l[ have proposed a model of gravity in which the mass of the 
graviton is taken into account. First attempts in this direction were proposed long time ago by Fierz and Pauli 0, 
who constructed a (ghost-free) linear theory of massive gravity. As later realized, the Fierz-Pauli theory is affected 
by some pathologies as, for example, it is in conflict with solar system tests 0]. A renew interest for this theory is 
arisen in the last years thanks to the Stueckelberg formalism introduced in Q (for a review, see [Hj]). 

The dRGT model relies on the idea to add higher order self-interaction graviton terms to the Einstein-Hilbert 
O h action in order to get rid of the Boudware-Deser instability 5] . An interesting consequence of the theory is that the 
| | presence of a tiny mass of the graviton gives rise to a term that play the role of cosmological constant. Therefore, 
(DO as a modified theory of gravity, it allows to account for the observed acceleration of the present Universe, without 
invoking exotic matter. In this respect, many efforts have been devoted to search for cosmological solutions of the 
field equations that describe the expansion of the Universe [l|, I(l-fl5| . For a Friedman- Robertson- Walker Universe 



ds 2 = a 2 (t) 



dt 2 - ^ - - r 2 (d6 2 + sin 2 Odcf 1 ) 
1 — Kr z 



(1) 



. where a(t) is the scale factor and K = ±1, is the spatial curvature, the modified Friedman equation (the — field 
00 " equation) reads 

o: 
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where the dot indicates the time derivative, k 2 = &Tt/m 2 P , mp = G- 1 ' 2 = 1.22 x 10 19 GeV is the Planck mass (in 
natural units), and ai, i = 0, 1, 2, 3 are constants 

a Q = c 4 - 4c 3 - 6 , at = 3C(3 + 3c 3 - c 4 ) , a 2 = 3C 2 (c 4 - 2c 3 - 1) , a 3 = C 3 (c 3 - c 4 ) , (3) 

where c 3 and c 4 are the parameters of the model, and C is a constant of integration. The i—j cosmological equations 
are not relevant for our aim, and therefore will be not reported here. The energy-momentum tensor of matter is 
chosen to be k 2 T^ — diag(p, —p, —p, —p) and satisfies the equation of continuity 

p + '3H(p + p)=0, (4) 

where H = a/a is the expansion rate of the Universe. Cosmological constraints on massive gravity theory have been 



recently studied in [17 



The aim of this paper is to study the massive gravity theory in relation to the formation of the light elements 
in the early Universe (big bang nucleosynthesis (BBN)). As well know, the big bang nucleosynthesis provides very 
stringent constraints on cosmological models, which must be satisfied in order to avoid conflicts with observations. By 
making use of the current observational data on primordial abundance of light elements, we infer constraints on the 
parameters characterizing the theory. As an application, we explore the consequences of the model in the framework 
of the thermal relic abundance. This analysis is particularly interesting because alternative cosmologies may modified, 
in principle, the thermal history of relic particles. This occurs during the pre big bang nucleosynthesis epoch, a period 
of the Universe evolution not directly constrained by cosmological observations. If the expansion rate of the Universe 
is enhanced as compared to the expansion rate of general relativity, then thermal relics decouple with larger relic 
abundance. The change in the Hubble rate may have therefore its imprint on the relic abundance of dark matter, 
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such as WIMPs, axions, heavy neutrinos. These studies have been motivated by the rising behavior of the positron 
fraction recently observed in the PAMELA experiment [l8l |. 

The paper is organized as follows. In Section II we derive the constraints provided by BBN that the parameters of 
the massive gravity theory must satisfy in order to be consistent with observational data. Section III is devoted to 
the analysis of dark matter relic abundance in the massive gravity cosmology. Conclusions are shortly discussed in 
Section IV. 



II. CONSTRAINTS FROM BIG BANG NUCLEOSYNTHESIS 



BBN provides, together with cosmic background radiation, a strong evidence that the early Universe was hot and 



dense. In the BBN epoch, the main interactions of particles are v e + n f>o+ e 
The weak interaction rate of particles in thermal equilibrium is given by [l9l . I20I 



A(T) ~ qT 5 



off 



(5) 



where A = A„ e+ „„ p+e - + A e 

++n<+p+u e + A.n-o-p+e-+& e and q — 9.6 x 10 46 eV 4 . Eq. ([5]) is obtained in the regime 
T> Q, where Q = m n — m p , and m„„ are the neutron and proton masses. To estimate the primordial mass fraction 
of 4 He, one defines the quantity [l9|, l20l | 



Y„ = X 



2x(t f ) 



<t f ) ' 



(6) 



where A 



-{U-tf)/r 



tf and t n are the time of the freeze-out of the weak interactions and of the nuclesynthesis, 



respectively, r ~ 887sec is the neutron mean life, and x(tf) — e~^ T ^f^ is the neutron to proton equilibrium ratio. 
The function A represents the fraction of neutrons that decay into protons in the time t £ [tf,t n ]. Deviations from 
Y p (generated by the variation of the freezing temperature Tf) are given by 



SY P = Y p 



Ye. 

2A 



In 
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(7) 



In the above equation we have set 5T{t n ) = because T n is fixed by the deuterium binding energy [21J. By making 
use of the current estimation on mass fraction Y p of baryon converted to 4 iJe during the big bang nucleosynthesis 

0E3 



Y p = 0.2476 . 



\sy p \ < nr 4 , 



one gets that the upper bound on ^j- is given by 



ST f 



< 4.7 x 10" 



(8) 



(9) 



Some comments are in order. Typically during the big bang nucleosynthesis the usual curvature term K/a 2 was 
unimportant, so that we shall set K = 0. Moreover, the Universe is radiation dominated with an energy density that 

scales with a 4 , 



p r ~ a , or in terms of the temperature, p, 



— - 9 * T A . The energy density p r is conveniently 
10~ 5 is the present value of fl r , 



30 

(0) (1 + z) 4 , where U° ] 



written in terms of the energy ratio fl r = y- = ^ 
p cr ~ 8.4/i 2 10 -47 GeV 4 is the critical density (0.5 < h < 0.8 is the normalized Hubble constant). 

Our aim now is to relate the variation 5Tf to the parameters © of the theory. We discuss the cases in which C ^ 
and C = 0. 



A. The case C / 

As shown in Eq. @, the ro 2 -corrections scales as a p , with p = 0,1,2,3. In the early Universe, the dominant 
term is given by a 3 -term. In the next, the latter will be treated a perturbation to the radiation energy density p r . 
Accordingly, the Hubble expansion rate H = a/ a given by Eq. ([2]) can be written as 




H = Jl + ^^H GR = H GR + I W ! - ^-4 I ) //„,,.. (10) 
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where Hqr = (n 2 pa 2 /3) 1 / 2 is the expansion rate of general relativity. 

The freeze-out temperature T = 7/(1 + ) follows by equating the Hubble expansion rate with weak interaction 
rate, H = A. One obtains 19] 



4 x AlATl = H GR {T f ) -> T/~0.6MeV, (11) 



and 



STf_ if m 2 a 3 \ m 2 a 3 

0.8\ 2 8.4/i 2 l(T 47 GeV' 



* LW ° (lo^Gev) UJ to ' (13) 

where we used the upper bound on the graviton mass m < 10~ 42 GeV and z ~ 10 9 . As expected, STf = in the case 
m = 0. 

Using the upper bound ([9]) one gets 

a 3 m 2 < 4.2 x 1CT 82 GeV 2 , (14) 



or 



2 



, ,'10- 42 GeV h V Per 



B. The case C = 



In this particular case, it follows that the contribution coming from massive gravity, i.e. m 2 ao (see Eq. ([2])), plays 
the role of cosmological constant. As discussed in Q, these terms may give rise to the late time acceleration of the 
Universe. 

The analysis goes along the previous one, with the replacement 1 + z — > (1 + z) 4 in (TT2|) . Therefore one obtains 



from which, using again ([9J , 



H 

Tf 



H, 



m 2 ao 
K 2 p 

m 2 ao 



GR 



10K 2 p cr n r 0) (1 + z) 



m 2 ao 
2k 2 P 



Hgr ■ 



(16) 
(17) 



a m z < 4.2 x 10" 55 GeV 



(18) 



Equations (| 14[) and (1181) represent the main results of the paper: they give the upper bound on the parameters 
characterizing the massive gravity theory that must be satisfied to have compatibility of the theory, and its predictions, 
with the present constraints dictated by the observational cosmology. Notice that these equations hold whatever is 
the mass of graviton. Given the latter, one can definitively fix the constants ao and a 3 . In what follows we shall refer 
to them for studying the PAMELA puzzles. 



III. CONCLUSIONS 



Massive gravity theory is an alternative covariant formulation of General Relativity. It relies on the fact that the 
graviton is massive and that negative energy states (ghosts) are absent. Interestingly, for an homogenous and isotropic 
Universe, the field equations turn out to be modified by a constant term that mimics the cosmological constant and 
by other terms that scale with the scale factor, mimicking dust, quintessence and stiff matter. 

Starting from the modified Friedman equation @ and using the recent bounds on the primordial light elements, 
we have derived the following upper bounds asm 2 < lCU 82 GeV 2 and aom 2 < 10 -55 GeV 2 . The former is obtained 
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assuming that in the early Universe the m 2 /a 3 -contribution is dominant with respect to the other corrections, while 
the letter is obtained setting C — (or assuming that the cosmological constant like term ~ a$m 2 is dominant). 

As final comment, we wish to discuss the interesting problem related to the recent results of the PAMELA ex- 
periment, i.e. the excess of positron events, which could represent a possible signal for dark matter through dark 
matter annihilation in our Galaxy. Theoretical results indicate that the PAMELA data can be understood if the 
annihilation cross sections are larger than those obtained in standard cosmology. A possible mechanism comes from 
alternative cosmologies, which give rise to an enhancement of the expansion rate of the Universe, hence large annihila- 
tion cross sections, being at the same time also compatible with other observations. To account for the enhancement 
of the expansion rate, it is usual to write 24-26] H(T) — A(T)Hgr(T), where A(T) is the enhancement function 

and H GR (T) = ^J^r Pr, with p r = ^^T 4 , the expansion rate of the Universe compute in the standard General 

Relativity. Moreover, conflicts with big bang nuclesynthesis predictions are avoided by working at temperatures of 
the Universe greater than the temperature at which the Hubble rate reenters the standard rate of general relativity 
(~ lMeV). In this regime, A{T) is conveniently parameterized as [24| 

A(T) = l + rj(^-J , (19) 

where Tp is a reference temperature at which the WIMPs dark matter freezes out in the standard cosmology (Tp ~ 
17.3GeV 24]). In general, Tp depends on the dark matter mass m x . rj and v are free parameters and characterize a 



specific cosmological model. The values of the parameter r\ required to explain the PAMELA data are 1 < r/ < 10 to 
which corresponds the WIMPs dark matter masses 10 2 GeV < m x < 10 3 GeV. For dark matter masses of the order of 
10 2 GeV, the parameter 77 can be also close to zero, i.e. m x ~ 10 2 GeV — y < r] -C 1. Assuming that the dominant term 

induced by massive gravity is a^/a 3 , and treating it as a perturbation of p r , one gets Eq. (|19p with 77 = a 3 16 ^l g "a"^ 

and v = -1 (T = 2.7K = 2.4 x lCT 13 GeV). Setting 7/ = 10~ x <C 1 it follows a 3 m 2 «6x lCT 72 - x GeV 2 , that is, 
comparing with the upper bound provided by big bang nucleosynthesis (|14j) . one needs x ~ 9 — 10, or r\ ~ io~( 9 ^ 10 ), 
in order that the massive gravity theory may explain the PAMELA experiment. The corresponding WIMPs mass is of 

the order 10 2 GeV. In the case in which C — 0, the amplification function (|T9|) is characterized by 77 = ao 16 ^l g "' T T F 

and v — —4. The value of 77 ~ 10~ 15 is required in order that results are consistent with (|T51) . The WIMPs mass still 
is ~ 10 2 GeV. 
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